Giant vortices in rotating electron droplets 
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We predict the formation of giant vortices in quasi-two-dimensional quantum dots at high mag- 
netic fields, i.e., in rapidly rotating electron droplets. Our numerical results of quantum dots confined 
by a flat, anharmonic potential show ground states where vortices are accumulated in the center 
of the dot, thereby leading to large cores in the electron and current densities. The phenomenon 
is analogous to what was recently found in rotating Bose-Einstein condensates. The giant-vortex 
states leave measurable signatures in the ground-state energetics. The conditions for the giant- 
vortex formation as well as the internal structure of the vortex cores are discussed. 
PACS: 05.30.Fk, 03.75.Lm, 73.21.La, 71.10.-w, 85.35.Be 
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I. INTRODUCTION 

In analogy to the well-known classical phenomenon, a 
vortex in a quantum system can be characterized by ro- 
tational flow of particles, forming a cavity at its center. 
It corresponds to a node in the wave function, associated 
with a phase change of integer multiples of 2n for each 
path enclosing it. The formation of vortices in quantum 
liquids has been observed, e.g., in superconductors jL He 
liquids^ and more recently, rotating Bose-Einstein con- 
densates (BECs) in atomic traps^ 

Recently, electronic structure calculations of quasi- 
two-dimensional (2D) quantum dots (QDs)^ have pre- 
dicted formation of vortices in rapidly rotating electron 
dropletsASii In this case an external magnetic field is 
used to set the system rotating. At sufficiently high fields 
the system tends to minimize the total energy by nucle- 
ating vortices, i.e., the quanta of the magnetic flux. This 
highly correlated quantum system of electrons and hole- 
like vortex quasi-particles 8 eventually shows formation of 
composite particles comprising an electron and a number 
of vortices^ The fractional quantum Hall effect of two- 
dimensional electron gas can be seen as bulk counterpart 
of this phenomenon^ The physics of rotating fermions 
in quantum dots shows remarkable similarities to that of 
rotating bosons in BECs>£ This can be seen as a conse- 
quence of the quantization of angular momentum in 2D 
systems. 

High angular momentum of a rotating group of parti- 
cles (bosons or fermions) is associated with the formation 
of multiple vortices which can form as a lattice or a clus- 
ter inside the system. Another possibility is that vortices 
are spatially concentrated into a small area at the center 
of rotation. This structure is characterized by a large 
localized core around which the phase shift is an integer 
multiple of 27T. We define this core as a giant vortex. 
Similar structures have been observed in experiments on 
rapidly rotating superfiow in BECs^ It should be noted 
that some authors reserve the term giant vortex for mul- 
tiply quantized vortices, 11,12 such as those found in the 



Laughlin states^ Our intention is, however, to analyze a 
phenomenon which is analogous to the formation of giant 
vortices in BECs. 

In this paper we report theoretical prediction that 
giant-vortex structures can be found in rotating quasi- 
2D fermion systems. These states are analogous to those 
found in BECs, and they emerge when anharmonicity 
(flatness) is added in the otherwise parabolic confining 
potential. Unlike in the bosonic case, the tendency of 
giant vortex formation is limited to systems with a low 
number of electrons (N < 20). As the second differ- 
ence, the cores of giant vortices in QDs do not merge 
completely. We show this effect using conditional wave 
functions to probe the internal structure of the many- 
electron state. We also suggest how the formation of 
giant vortices could be observed in experiments. 



II. MODEL 

Our model system is a QD with electrons being re- 
stricted in quasi-2D plane and interacting via their mu- 
tual Coulomb repulsion. The Hamiltonian is then 
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where N is the number of electrons, V c is the external 
confining potential on the plane, and A is the vector po- 
tential of the homogeneous magnetic field oriented per- 
pendicular to the QD plane. We imitate the conditions of 
a real semiconductor heterostructure by using the effec- 
tive material parameters for GaAs, i.e., the effective mass 
to* = 0.067 m e and the dielectric constant e = 12.4 eo- 

We solve the many-particle Schrodinger equation for 
the Hamiltonian Q using two numerical approaches, 
namely the configuration interaction (CI) method which 
numerically diagonalizes the Hamiltonian, and the spin- 
density-functional theory (SDFT) in a real-space formu- 
lationiii We point out that our previous studie oV 5 show 



that SDFT provides a proper description of the high-field 
solutions even if the total energies given by the computa- 
tionally more demanding current-spin-density-functional 
theory are slightly closer to the quantum Monte Carlo 
results— 

According to theoretical studies of rotating BECs, 
giant-vortex structures may emerge in bosonic systems 
when the confining potential V c has a flat, non-parabolic 
shape— Therefore, within the CI method we apply 
a confinement following the notation of Jackson and co- 
workers^ 
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Here Hloq is the confinement strength fixed to 5 meV un- 
less stated otherwise, a = ^h/m*uo ~ 15.1 nm is the 
oscillator length, and A is a positive dimensionless con- 
stant determining the strength of the quartic contribu- 
tion, i.e., the flatness of the confinement. Within the 
SDFT approach we apply also a circular hard-wall po- 
tential well which can be considered as the extreme limit 
of the flatness of the potential. It is defined simply as 
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where R is the radius of the well. It should be noted that 
as a reasonable approximation, the confining potential in 
semiconductor QDs is usually assumed to be parabolic, 
although detailed modeling of realistic QD systems has 
pointed out the importance of anharmonicity in the po- 
tential^ In addition, steep edges have been found to 
explain observed features in far-infrared absorption spec- 
tra— As we shall see below, even small deviations from 
the harmonicity may lead to emergence of new types of 
ground states. 



III. PHASE DIAGRAM 

First we visualize the fundamental difference between 
the high-field solutions corresponding to a vortex cluster— 
and a giant vortex in a QD. Figure ^ shows the SDFT 
results for electron density of a vortex-cluster state (a) 
and a giant-vortex state (b) in a 10-electron QD defined 
by a parabolic (A = 0) and circular hard-wall confining 
potentials, respectively. In the parabolic case the vortic- 
ity is spread throughout the system which gives rise to a 
cluster of seven hole- like quasi-particlesi^ If the confine- 
ment is set flat, however, we find a large core in electron 
and current densities at the center. In this case all phase 
singularities are absorbed into the core region to form 
a single seven-fold giant vortex. The qualitative origin 
of the phenomenon is in the fact that increasing anhar- 
monicity (flatness) in the external confinement leads to a 
relative increase in the energetic advantage of the outer 
orbitals compared to the inner ones. In consequence, 
there is rearrangement in the electron occupations such 




FIG. 1: (Color online) Electron densities (SDFT) of N = 10 
quantum dots containing seven flux quanta in both but having 
different distributions of vorticity. The state in (a) is calcu- 
lated using parabolic external potential (hui = 4 meV, B — 13 
T) and shows a cluster of seven vortex holes. The state in (b) 
is calculated in hard-wall potential well (R — 100 nm, B = 25 
T) and shows a giant core (eye) of a vortex comprising seven 
flux quanta close to the center of the rotation. 



that the inner orbitals are left empty (or they have a 
negligible weight). Correspondingly, the vortices (holes 
in the electron density) are accumulated at the center. 

To gain more physical insight into giant vortices, we 
present in the following a detailed analysis of the vortex 
states in a six-electron QD. In a parabolic potential the 
calculations have shown that the ground states in the 
fractional quantum Hall regime (filling factor v < 1) oc- 
cur only at certain magic values for the angular momen- 
tum L (Refill). For N — 6 the v = 1 state corresponds 
to L = 15 maximum-density droplet (MDD) and the sub- 
sequent magic angular momenta are L = 21, 25, 30, 35,... 
These states correspond to ground states with an increas- 
ing number of off-electron vortices forming clusters in- 
side the electron droplet £ When the confining potential 
is made flatter by including the quartic contribution of 
Eq. (J2J , new ground states emerge corresponding to giant 
vortices. The CI calculations show that when A > 0.06 a 
double vortex core appears between L — 25 and L = 30 
states and it has an angular momentum L = 27. Triple 
vortex core emerges for A > 0.16 and it has L = 33. This 
is illustrated in Fig.^a) which shows the phase diagram 
of the CI solutions for < A < 0.2. It should be noted 
that the minimum value of the quartic multiplier A re- 
quired for giant vortices is relatively small (A rj 0.06), so 
that the results may apply to a large class of real quantum 
dot devices. Fig. Ufa) shows also that the giant-vortex 
states vanish rapidly at high fields unless A is consider- 
ably increased. This is due to the increasing magnetic 
confinement which has a parabolic form. 

Fig- Etb) shows the magnetization M — —dE to t/dB 
of a QD with A = 0.18 associated with the dashed line 
in Fig. [21 a). The steps in M correspond to changes in 
the many-electron ground state either due to an increase 
in the vortex number— or due to a rearrangement of the 
vortices (and electrons) in the system. The oscillative 
behavior of the magnetization might be observable in, 
e.g., magnetization measurements for large ensembles of 
QDs. Such experiments can been done using sensitive mi- 
cromechanical magnetometers— Our data suggest, how- 
ever, that the ground-state transitions associated with 
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FIG. 2: (Color online) (a) Phase diagram of the vortex so- 
lutions in a six-electron quantum dot defined by a parabolic- 
plus-quartic confining potential [Eq. J2J]. The results have 
been calculated with the CI method. The Roman and Arabic 
numerals mark the number of multiple and clustered vortices 
in the ground states, respectively. The total angular mo- 
menta L are marked by circled numbers, (b) Magnetization 
of a quantum dot with A = 0.18 corresponding to the dashed 
line in (a). 



the rearrangement of vortices and electrons usually cause 
relatively small changes in the magnetization. Hence, the 
observation of such transitions would require a high de- 
gree of accuracy in the experiments. 



IV. INTERNAL STRUCTURE OF GIANT 
VORTICES 

In order to visualize the internal structure of the giant- 
vortex states, we compute the conditional wave func- 
tions^ from the numerically exact many-body ground 
states. Within the SDFT instead, we use the auxil- 
iary single-determinant wave function constructed from 
Kohn-Sham states. 25 ' 26 Figure |31 shows the total elec- 
tron and current densities (left panel) and the conditional 
wave functions (right panel) of the double (a) and triple 
(b) vortex states at L = 27 and L = 33 (CI) and a giant 
eight-fold vortex state (c) in a circular well with L = 63 
(SDFT). Since the system is totally polarized there is a 
Pauli vortex on each electron position as mandated by 
the exclusion principle. The exchange-hole around each 
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FIG. 3: (Color online) Left panel: Total electron densities 
(grayscale) and current densities (arrows) of double-vortex 
(II), triple- vortex (III), and eight- fold giant- vortex (VIII) 
states in a six-electron quantum dot. The states II and III 
have been calculated with the CI in a quadratic-plus-quartic 
potential as in Fig. [5] The state VIII has been calculated with 
the SDFT in a circular hard-wall well of radius R = 70 nm. 
Right panel: Corresponding conditional densities (contours) 
and the phases (grayscale) of the conditional wave function. 
The phase of the SDFT result (VIII) is calculated from the 
single-determinant auxiliary wave function constructed from 
the Kohn-Sham states. At the lines where the shadowing 
changes from darkest gray to white the phase jumps from 
— 7T to 7T. The crosses (blue) and circles (red) mark the most 
probable electron positions and vortex positions, respectively. 
The Pauli vortices at the electron positions are not marked 
by circles for clarity. The probing electron has its maximum 
probability downright. 



electron can be associated with these vortices. In addi- 
tion, there are off-electron vortices close to the center of 
rotation. They lead to holes in the electron density with 
rotational current around them. The size of the vortex 
hole at the center of the QD increases with the number 
of vortices as shown in the left panel of Fig. [3] 

Despite the formation of single cores in the electron 
and current densities, the conditional wave functions 
[right parts of Figs. EJa) and (b)] show finite separation 
between the vortices at the center. The separation is 
due to quantum fluctuations appearing as effectively re- 
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pulsive vortex-vortex interactions. This effect can also 
be interpreted as zero-point motion of vortices^ which 
leads to non-zero electron density at at the center of the 
giant vortex. In the bosonic case the vortex- vortex inter- 
action energy as a function of vortex separation has been 
approximated to be of logarithmic form which prevents 
complete merging of the vortex cores m& 

In the SDFT result for the eight-fold giant vortex 
[Fig.JSfc)], the density is zero at the center and very small 
up to the QD radius where the electrons are strongly lo- 
calized. The vortices at the center are then packed close 
to each other but still, however, showing a small-scale 
structure. The stronger vortex localization in the SDFT 
can be understood from the fact that the SDFT gener- 
ally can incorporate less correlation effects than exact 
many-body methods. 

Qualitatively similar internal structures have been 
found for giant vortices in rotating BECsi-^ In that sys- 
tem, however, a local potential maximum at the center 
attracts the phase singularities into a single point in space 
so that a multiply quantized vortex is formed^*2£ This 
is not the case in QDs, where the confining potentials 
increase monotonically as a function of the radius of the 
system. In QDs a pinning potential, expectedly leading 
to perfect multiple phase quantization like in bosonic sys- 
tems, could be realizable in form of a negatively charged 
impurity at the center of the sample. 



V. DEPENDENCE ON THE ELECTRON 
NUMBER 

When N is increased the maximum number of flux 
quanta that can be packed to the center of a giant 
vortex decreases rapidly. According to our SDFT cal- 
culations for hard-wall QDs [Eq. @] of realistic radii 
(R = 50... 150 nm), giant vortices consisting of more 
than ten flux quanta can be found when the electron 
number is increased up to N = 10. When N = 12, 
however, only three flux quanta can be packed at the 
center. For N > 18 giant vortices disappear completely 
for realistic well sizes, and they are replaced by vortex 
clusters inside the electron droplet. We visualize this ef- 
fect in Fig. which shows the total angular momenta 
of the SDFT solutions for N = 12 (a) and N = 20 (b) 
hard-wall QDs as a function of the magnetic field. Red 
and blue dashed lines mark the L plateaus for multiple 
vortices (Roman numerals) and vortex clusters (Arabic 
numerals), respectively, and the black dashed line corre- 
sponds to the single- vortex solution. The results between 
the plateaus are a consequence of mixing of the consecu- 
tive exact ground states^ In the range of the MDD and 
the highest multiple vortex solution, the differences in L 
between the plateaus are AL = N. This regularity fol- 
lows from the simple increase in the hole size discussed 
above. When the vortices become clustered instead, we 
find smaller steps in L resulting from the complex rear- 
rangements of vortex clusters in the QD. 



(a) 

v '67 
56 
-q 45 

Q 

J 5 36 
i 

=± 24 
12 



6_. 
Ill 



AL = N-1 



| N-1 
J N-3 



II 



1 = 1 

MDDy 



10 



12 



14 16 
B [T] 



18 



20 




B [T] 

FIG. 4: (Color online) SDFT results for the total angular 
momenta of N = 12 (a) and N — 20 (b) hard-wall quantum 
dots. Red and blue dashed lines mark the L plateaus for mul- 
tiple vortices (Roman numerals) and vortex clusters (Arabic 
numerals) , respectively. The black dashed line corresponds to 
the single- vortex solution. 



The disappearance of giant vortices as a function of N 
is related to the formation of shell structure in the QD, 
i.e., it becomes energetically favorable for the electrons 
to occupy the inner orbitals of the system. This behav- 
ior is in striking contrast with what is found in bosonic 
systems where giant vortices can be found in BECs with 
a very large number of particles. On the other hand, our 
SDFT results are consistent with the exact diagonaliza- 
tion results which show that even single vortices tend not 
to localize at the center of a parabolic QD when N > 12 
(Ref.lH). 



VI. SUMMARY 

To summarize, on the grounds of our numerical anal- 
ysis we predict that giant vortices emerge in quantum 
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dots in high magnetic fields, or more generally, in rapidly 
rotating electron droplets. These ground states are anal- 
ogous to those found in rotating Bose-Einstein conden- 
sates. In quantum dots only a slight anharmonicity (flat- 
ness) in the external confining potential is required for 
the formation of giant vortices. By using conditional 
wave functions we can analyze the internal structure of 
the vortex cores. Formation of giant vortices leave mea- 
surable signatures in the ground state electron and cur- 
rent densities as well as in the magnetization of quantum 
dots. These signatures can, in principle, be detected in 
state-of-the-art magnetization or electron transport mea- 
surements or via direct imaging of the charge density in 
quantum dots. As in the bosonic case, the conditions for 
the emergence of giant- vortex states are determined by 
the interplay between the interactions and the form of 
the confining potential. We find that in quantum dots 
the giant-vortex solutions are generally limited to sys- 
tems consisting of less than 20 electrons. However, our 



numerical analysis has been done assuming the conven- 
tional electron-electron interaction of a Coulombic form. 
Other types of interparticle potentials may lead to quan- 
titatively different results, but we leave this topic for fu- 
ture research. 
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